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The quality that can be defined is not the Absolute Quality 
The names that can be given to it are not Absolute names 
It is the origin of heaven and earth 
When named it is the mother of all things ... 

Tao Te Ihing 
of 

Lao Tzu 


GHAPTER I 
IRTROLUOTIOR-. 

1 .1 POWER SYSTM STABILITY 

The term 'Stability* when used in the conteirf; of a 
power system is that property which ensures that the system 
will remain in operating equilibrium through the normal and 
abnormal operating conditions [13. As the structure of power 
systems grows larger and more complex, the stability studies 
gain paramount importance as far as the safety of the total 
system is concerned. With the ever increasing demand for 
electrical energy and dependence on an uninterrupted supply, 
the associated requirement of high reliability dictates that 
power ^sterns be designed to maintain stability under most 
disturbances, consistent with economy. Furthermore, the 
degree of stability of a power system, is an important 
consideration in the planning and designing of new facilities. 

Investigation of the disturbed conditions in a power 
system following a fault is carried out on a mathematical ‘ 
model of the ^st«n. It is customary to classify stability 
[l,2,5l, for the purpose of analysis as ‘Steady state stability’ 
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and ’Transient stability*. The occurrences disturbing the 
stable operation of a power system can be small and gradual 
changes, as for instance controlled changes of power input 
and output, or large changes in operating conditions such 
as unco itrollable loss of li.rge amounts of generation or 
load, or the occurrence of heavy tranani ssion line faults. 

The terms gradual and sudden should be interpreted here in 
relation to the major time constants of the system or part 
of the system under study. Small and gradual changes 
constitute problems in steady state stability while large 
and sudden changes lead to problems in transient state 
stability, 

1.2 TRAWSIMT STABIIITI 

Transient stability as referred to a single synchro- 
nous machine connected to an infinite bus may be defined [4] 
as follows. 

nE^;IiriTIO]T: ’If a synchronous machine operating in steady 

state torque equilibrium issubjected to a disturbance of 
any kind which results in speed deviations (oscillations) of 
the machine rotor from the (synchronously rotating) reference 
axis, then the machine is called transiently stable if the 
rotor swings die out and the rotor reaches a new stable 
equilibrium position’. 

Transient stability for a multimachine system may 
be defined [4 ] as follows. 
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DEPIITITIOir: 'If individual machines in a multimachine systoE 

are operating in a steady state torque equilibrium and a dis- 
turbance of any kind is imposed on the system, then the 
system is called transiently stable if each machine oscillates 
around and ultimately comes to rest at a new stable torque 
equilibrium point’ , 

The large perturbation which creates the transient 
stability problem, may be sudden increase in load, or a 
sudden increase in reactance of the circuit, caused for 
instance by a line outage. But the most severe perturbation 
to which a power system can be subjected, is a short circuit, 
generally termed as a ’fault', and a symmetrical three phase 
short circuit is the most severe of all faults. In designing 
and planning a system all possible contingencies including 
the worst conditions have to be considered. Hence in 
stability studies a three phase short circuit is generally 
considered [23. 

In faulted systems, it is essential to determine when 
the fault has to be cleared without the rotating machines 
losing synchronism. The maximum time through which the 
fault can be left on the system so taat the rotating machines 
will return to their normal conditions after the fault is 
cleared (by opening one or more circuit breakers) is called 
the ’ critical clearing time’ . Determination of critical 
clearing time constitutes an important aspect of transient 



stability study since it is this infomation 'Which enables 
the setting of relays, -which in turn control the tripping of 
the circuit breakers. 

1 . 5 GOWBJTIOFAI METHODS POR TRAHSIEIIT STABILITT ANALYSIS 

The analysis of transient stability is formulated as 
follows: Given a system does there exist an equilibrium 
state of the system after the dis-turbance is cleared? If 
yes, what is the maximum time that the dis-turbance may be 
allowed to remain without the ^stem losing synchronism? 

The analysis consists of two main steps. 

STEP 1: The study of the evolution of the system from the 
occurrence of the dis-turbance to the time of 
clearing. 

STEP 2: The study of the evolution of the system beyond 
the time of clearing. 

The intervals of time corresponding to these -two steps are 
called the ’faulted state' and the 'postfault state* respec- 
tively, Each of these steps involves th© integration of a 
large number of nonlinear differential equations. 

The classical method of stabili-ty investigation 
consists in assuming a clearing time ' t^* arbitrarily and 
solving the differential equations of the system in its 
faulted state upto t^, and in the postfault state after t^ 
to obtain the swing curves. The initial eonditions are 
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olstained from a prefault load flow data. The swing curves 
which are a plot of the rotor angles versus time are examined 
for the stability of the postfault system. This is usually 
done by observing whether the rotor angle differences with 
respect to any one machine taken as a reference tend to 
remain constant during the first swing of the rotor an^es. 

In such a case we conclude that the system is stable. If 
they tend to diverge, then it is considered to be unstable. 
Generally it is assumed that if the system is stable during 
the first swing, it is also stable thereafter. This process 
is repeated for another clearing time. After a certain 
number of such repeat procedures, the clearing time at which 
the transition from the stable to the unstable state of the 
system takes place is established. This gives the critical 

clearing tind This pracodurd of d'otormination of "ttie 

critical clearing"' time is not efficient computationally when 
one requires a quick picture of the stability with respect 
to various possible perturbations. It is because of this 
reason, that special onphasis has been laid on examining 
alternative approaches such as direct analytical methods, 
in the last decade. 

Stability investigations by direct methods in 
contrast to the classical methods are carried out by 
making partial or no use of the system differential 
equations. The equal area criterion [ 2 ] is a graphical 


II 



6 


procedure for single and two machine systems and does not 
make use of the differential equations as such. Analytical 
methods such as the phase plane technique [5] were employed 
for single and two machine systems. The energy integral 
criterion [6 J and the direct method of lyapunov [7-23 ] for 
multimachine systems make use of the system differential 
equations partially. The Lyapunov’s direct method has been 
a topic of active research in recent years. This method 
essentially consists in establishing a region of stability 
around the postfault equilibrium state via the construction 
of a Lyapunov function. The system is considered stable if 
the solution trajectories of the faulted and the postfault 
states lie entirely in this region [ S ]. 

1 . 4 APPLIOATIOTir OP LYAPUN-QY’ S LIREOT METHOD TO THE 
TRAHSIEHT STABILITY PROBLEM OP POWER SYSTMS 

The crux of the problem, in the application of 
Lyapunov’s direct method to the transient stability analysis 
of multimachine power systems, revolves around the replace- 
ment of the step 2 detailed in Section 1.5* by a stability 
criteritLi in the form of a Lyapunov function around the 
postfault equilibrium point, Aylett [6 jin 1953, was the 
first to suggest such a criterion based on energy cwnsidera- 
tions, but the theory was not sound enou^ for practical 
implementation. Pioneering work on the application of 
Lyapunov’s direct method to power system transient stabi- 
lity was done by &1 ess [7 ] and H-Abiad and Hagappan [sj. 
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The approach consists of writing the system differential 
equations in the postfault state (after the final switching) 
in the form 

X = f(x) (1.1) 

with the postfault equilibrium state at the origin 0, and 
seeking a suitable Lyapunov function ’V'(x) which, along with 
its time-derivative 'V(x) has the required sign-definite 
properties. The stability region around the postfault 
equilibrium state 0 is specified by an inequality of the 
type V(x)<0, where 0 is a constant to be determined. This 
constant is usually obtained by evaluating the function 
V(x) at the unstable equilibrium point closest to the 
postfault stable equilibrium point. Wow the critical 
clearing time can be calculated by numerical integration 
of the faulted system differential equations, substituting 
at each step the value of x into V and checking whether 
V < 0, When V reaches the value 0 (i.e., the faulted system 
state trajectory reaches the boundary V = 0), the time at 
that instant is the estimated critical clearing time ' • 
Thus an explicit integration of the differential equations 
is performed only once, resulting in a narked reduction in 
the computation time, especially in the analysis of large 
scale power systems for various contingencies. Another 
advantage of this method lies in investigating, the effects 
of parameter variations such as damping [9,26], power input 
[27 l etc, r 
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It should however be noted that Lyapunov theorems 
provide only sufficient conditions, and not ’necessaiy and 
sufficient’ conditions, for stability. Thus, the estimation 
of stability regions and critical clearing time based on 
Lyapuno ’s direct method may turn out to be conservative, 
but it can be used .very effectively as a preliminaj;^ 
screening procedure, and reduce the number of studies to 
be carried out in detail. 

1 .5 THE LARGE SCALE POWER SYSTEM STABILITY PROBLIM 

The previous section explains the technique of 
applying Lyapunov methods. Although, concentrated efforts 
have been made towards the improvement of the Lyapunov 
functions [12,17], the application of the theory to realistic 
systems has been very slow. This is due to the fact that 
there are a. few drawbacks, apart from the inherent 
conservativeness of the method. They are (i) The applica- 
tion of the theory to large power systems poses heavy 
computational problems in the determination of staMlity 
region^ based on computing T(x) at the nearest unstable 
equilibrium point (u.e.p.) [2S]. It has been established 

that for a n-machine system, the number of unstable equili- 

i 

brium points is 2’^“ -1, and consequently very large when ’n* 

is large. This problem has received considerable attention 
re cently [ 21 , 22] , but requires further refinement to obtain 
more acoirate stability regions, (ii) The selection of the 
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models for the system, is strongly constrained by the 
possibility of a finding a suitable Lyapunov function [28], 

It has not been possible to include refinements like accurate 
models of synchronous machines, saturation effects etc., into 
the Lyapanov functions for miltimachine systems. Most of the 
work is Confined to classical models with other effects 
included partially, (iii) As size of the system grows the 
overall computational burden in using Lyapunov functions 
also increases. 

The above considerations therefore call for new 
theoretic approaches to the problem. Two such promising 
approaches are: 

(i) Dynamic equivalencing [29-34] 

(ii) Use of vector Lyapunov functions [36-59], 

The first approach consists in obtaining a simplified model 
after coherency analysis and dynamic equivalencing. Lyapunov 
functions are then constructed for this reduced order model. 

The second method uses a theoretical approach which 
was first proposed by Bellman [3?] and Bailey [361. They 
demonstrated its usefulness for the stability analysis of a 
complex composite system. Most complex systems are made up 
of a number of low order interacting systems. These low 
order systems in hierarchial theory [38] are called 
’ subsystems’ . These subsystems are arrived at on the basis 
of decomposition of the large system. 
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Today’s power systems being large and complex become natural 
candidates for the application of the decomposition - aggrega- 
tion technique. The teclmique consists in decomposing the 
composite ^ystan into low order subsystems. Suitabl e Lyapunov 
functions are constructed in the usual manrer for these 
systems, satisfying the required sign definite properties* 

On a higher hierarchial level, the concept of a vector 
Lyapunov function is introduced in which these scalar 
functions form the elements of the vector Lyapunov func- 
tion defined for the composite system. Using the proper- 
ties o f in teractions among the subsystems, conditions and 
region of stability for the overall system are derived. 
Initially this technique was restricted to a particular 
class of interconnections viz. linear interconnections. 
Nonlinear interconnections were proposed by Piontkovskii 
et al. [39]. Recent works of Siljak [40-43], Grujic and 
Siljak [ 44 - 45 ], Thompson [7^, Araki and Zondo [47], 

Michel [48-50], Pai and Narayana [51I Jocic, Ribbens 
Pavella and Siljak [52], and Jocic and Siljak [53] give 
an explicit exposition of these concepts, and provide an 
impetus to apply the results to large scale power systems. 

1.6 TRANSFER OONDUOTANOES IN MULT IMA CHINE POWER SYSTEMS 

In the application of Lyapunov’s method to realistic 
systems there arises a difficulty of taking into account the 
effect of the transfer conductances contained in ihe system. 
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In early investigations, transfer conductances -were 
entirely neglected or only partially taken into account [8,12]. 
Althou^ the conductance components take small values as 
compared with the susceptance components, it still remains 
an open problem whether or not we can neglect the conduc- 
tance components from the standpoint of the transient stabi- 
lity. Por certain cases, as shown by Ribbens-Pavella [14], 
the neglect of the transfer conductances may be justified. 
However, that may not be the case in general. Especially, 
they cannot be neglected if the system is heavily loaded. 

It has been pointed out by Uemara et al. [60 ] that there 
exists a danger in judging a practically unstable system 
to be stable if we use a Lyapunov function which neglects 
the transfer conductances. 

In recent works [61-62] methods have been proposed 
of including transfer conductances for transient stability 
analysis. 

1.7 OBJECTIVES AHD OUTLIITES OF THE THESIS 
The objectives of the thesis are: 

i) Application o f Lyapunov methods to large scale power 
systems using dynamic equivalencing. 

ii) Investigation of new decomposition structures for 
vector Lyapunov function application to power systems and 
validating on a raultimachine power system. 

iii) Examining effects of transfer conductances on 
the critical clearing time. 
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In Chapter II, Lyapunov's method for a realistic power 
system consisting of 13 machines, 71 "buses and 94 lines is 
applied for the full system as well as for the reduced system 
obtained by dynamic equivalencing procedure. For the latter, 
the system is divided into the study system and the external 
system. Coherency analysis is then done for generators in 
the external system . The coherent groups are egair^ised using 
the method proposed by Podmore [29 1 Lyapunov method is 
then applied to the reduced system. Critical clearing times 
are computed for various fault locations in the study system 
and the results obtained are compared with those using 
Lyapunov method for the full system as well as with those ^ 
obtained by a base case transient stability analysis. 

Chapter III presents a new decompesition technique 
based on the' centre of angle' principle proposed by 
Stanton [54] and Tavora and Smith [55], for the application 
of vector Lyapunov methods to large scale power systems, A 
vector Lyapunov function for the composite system is then 
found, and an overall stability region is constructed using 
the separate results of "Waiter and MoClamroch [59], 
Weissenberger [56], Grujic and Siljat [45], Araki [57] and 
Michel [48]. A four machine power system considered earlier 
by El-Abiad and Fagappan [8] is worked out using this method 
and the results are quite en couraging. 

In Chapter I"7 Lyapunov functions containing partially 
the effects of the transfer conductances, as proposed by 
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Pai and Varwandkar [ 62 ] are used in stability analysis of 
the UPSEB system, and the critical fault clearing times are 
calculated. Oomparison is made with the critical clearing 
times obtained without the transfer conductances. 

In the concluding chapter, the results of the thesis 
are reviewed and problem for future research outlined. 



Your knowledge is your friend in distant lands. 

Your wife is your friend at home. 

To the side, the right medicine is the friend. 

Dharma is a friend even "beyond the graves. 

»TEmS» 

CHAPTER II 

DYNAMIC EQUIYAIENCING AND APPLICATION OP LYAPUNOV’S 
METHOD TO A PRACTICAL POWER SYSTEM 

2.1 INTRODUCTION 

During the last decade, several research workers have 
applied Lyapunov’ s direct method for the analysis of the 
stability of power systems following a sudden distur"bance. 

It still continues to "be a fertile area for research for 
power engineers all over the world "because of the yet unsolved 
pro"blems seeking solutions. One of the pro"blems is the 
dimensionality of the large scale systems. In spite of 
improved Lyapunov functions "being availa"ble, computing 
effort is still aonsidera"ble due to the need for computing 
efficiently the stability regions. These regions are computed 
by evaluating V(x) at the nearest u.e.p. One of the ways to 
reduce computational time in evalua ting u . e.p. is through 
dynamic equivalencing procedures and model reduction tech- 
niques. 

In this chapter Lyapunov’ s direct method is applied 
to the UPSEB system consisting of 13 machines, 71 buses and 



94 lines using dynamic equivalencing procedure. Coherency 
analysis is first done for the external ^stem, and the 
group of coherent generators is found. These generators are 
then equivalenced using a method proposed by Podmore [29-31 
Lyapunov’s method is then anplied to the r'^duced order model, 
and critical clearing times are evaluated for a' set of 
representative fault locations in the study system. Critical 
clearing times for the same set of faults are computed for 
the original system.? using lyaxunov's method. A comparison 
of the two results, validates the dynamic equivalent, and also 
helps to prove the feasibility of applying lyapunov’ s method 
to reduced order models. There is a considerable saving in 
computer time through this procedure, 

2.2 MATHEMATICAL MODEL OP THE POWER SYSTEM 

A number of simplifying assumptions are made in 
arriving at a mathematical model for whi ch satisfactory 
Lyapunov functions can be constructed. These are: 

(i) The transmission network is considered to be in 
steady state which results in an amenable analytical expres- 
sion for the electrical power P^^ delivered by eadi generator. 

(ii) The voltage behind the transient reactance is assumed 
constant. In other words it is assumed that the flux linkages 
of various madiines are constant. 

(iii) Damping is assumed to be directly proportional to 
the slip velocity and is thus mainly due to the mecSianical 
friction and asynchronous torques. 
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(iv) The mechanical power input to each machine during 
disturbance is assximed constant and equal to the input 
prior to occurrence of the fault. This assumption is 
justified because the governor time constants are much larger 
than the duration of the transient swings. 

(v) loads are represented as constant impedances which 
results in the overall -'nodal admittance matrix at the internal 
nodes of the generator being constant. 

(vi) Transfer conductances in the admittance matrix are 
usually neglected, 

(vii) Saturation is neglected. 

(viii) Only round rotor machines are considered thus 
neglecting the effect of saliency. 


Under the above assumptions the equation governing 
the rotor motion of each machine is given by 


M d 6i 
^ dt2 


d6. 


+ 1 • 


i dt 



i=1 , 2 , . . . ,n 


( 2 . 1 ) 


where 

- inertia constant of th' ith machine 

- damping coefficient of the ith machine 

6^ - rotor angle measured with respect to a synchronously 

rotating reference frame 

P . - mechanical power input to the ith machine 
mi 

P . - electrical power delivered by the ith machine. 
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The electrical power delivered by the ith machine is 
given by 


ei 


4 ®ii 


n 


3=1 

3?^i 


S. Y. . 
3 13 


0os((6^ - 5 .) - 9. ) (2.2) 


3 ' '13 

i=1 , 2, . . . ,n 


wh er e 

s= - internal voltage 

^ij “ transfer admittance between the ith 

and the jth machines 

—1 

9j_lj = "tsin ) - the admittance angle between the ith 

and the jth machines. 


Substitution of (2,2) into (2.1) yields 


d^6 


d6 . 


1 T\ i -ri •m2 


+ D. At = 

1 JJ.2 1 dt mi 1 11 

dt 


- A- ®3 ^13 


3=1 

35^1 


i=1,2,...,n (2.3) 

Use of the appropriate values of 6-- ^ in (2.3) 

render them valid either for the faulted or the postfault 
state. The postfault steady stated ? is given by the 
solution of (2.3) with ^ » 0,*^ pertaining to the postfault 
network conditions and 6^ = 0, 6^ = 0. This implies that 
the angles 6? (i = 1,2,...,n) constitute a solution of the 
equations 


^mi"®i ^ii 


!^i 


E. E. Y. . los( 6.. - 0. .) i=1,2, 

1 3 13 13 13 ’ 


,n 


(2.4) 
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Since satisfies Eq.(2.4) the swing equations (2.3) can he 


expressed as 

d^6 d6. 

M. - 5 + D. , . ^ 
dt^ 1 dt 


n 


= - C % E. Y,. 0o.(6. .-e.j)-0os(a°,-6^^) 

i=1,2,...,n (2.5) 


Expanding the right hand side trigonometrically and using the 
relation Gos 9^^. and Sin ^6 ge'fe 


d^6. d6 

M.. 2 + 

i^t^ 1 dt 


n 

c 

D=1 

dA 


E. 

1 


E. 


G - ^ ( Go s 

J- tj 




+ B.^(Sin d^j-Sin 6°^) , i=1,2, . 


»n 

(2.6) 


The equation (2,6) can he expressed in the following 
state variable form 


X = Ax - B f(cr) 
G X 


(2.7) 


The formation and structure of matrices A, B and G for any 
given system is discussed i i references [21 ,32 3. 

2.3 STEPS I^TVOIVEI) HT APPLYING LYAPUNOV'S METHOD TO A 
POWER SYS'TEM 

The procedure for applying Lyapunov's method to a 
power system and determining the critical clearing time, 
consists of six main steps. 
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1. Load flow for the prefault system, 

2. Determination of driving point and transfer admittances 
between the internal buses of the machines for the 
faulted and postfault system, 

3. Determination of the postfault stable equilibrium point 

-if 

and testing its stability. 

4. Oons traction of a Lyapunov function (V..fun ction) about 
the postfault stable equilibrium point of the system, 

5. Determination of the boundary V 0 of the stability 
region, where 0 is the value of the V-funotion at the 
unstable equilibrium point closest to the stable 
equilibrium point for the postfault system. 

6. Integration of the faulted systan equations, until 
tragectoiy reaches the boundary. This gives the critical 
clearing time. 

2.4 ■ THE STABLE. EQUILIBRIUM POINT 

This is a very important step in the formulation of 
the stability criterion. The postfault systan differential 
equations are given by Eq,(2,1), Equations (2,1) can be 
cast in the following state-- variable form. 

X = E^Cx) (2.8) 

The equilibrium states of the postfault system can be found 
by solving the following nonlinear algebraic equation 

E|(x) = 0 (2.9) 
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The stable solution of (2.9) is given by minimizing 
vicinity of the prefault stable solution by 
any ■well known method like the Pletcher-Powell or oonjugate 
gradient method. The minimizatio’i problaa beaDmes time 
consuming and complicated (in the sense of convergence) as 
the size of the system (i.e. number of state variables) 
increases. To eliminate the unwl^dy minimization procedure, 
a new special load flow techni^lue proposed in ref, [24] has 
been used to compute the postfault stable equilibrium point. 


Having reduced the system to generator internal nodes, 
each node except one which is the slack node is treated as a 
’voltage controlled node’ (P-V node). The node which 
corresponds to the machine having maximum generation at 
the prefault stage is generally chosen as the slack node. 

The solution of the load -flow problem is initiated by 
assuming voltage angles at each node equal to the prefault 
steady state angles, and voltage magni'tudes equal to the 
scheduled voltages. Therefore voltage at the pth node will 


be 

+ 3 Sin S^) ; p=1,2,...,n (2.10) 


where pj is. the 

com pi ex power flowing 


scheduled voltage at the pth node 
into the pth node is given by 


4 


The 


P - 3Q = V* X]; V_ 
P P P ^ PI '1 


( 2 . 11 ) 


It is to be noted in equation (2.11) that P.^^ is equal to P 


mp 
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i.e, power delivered by the machine connected to the pth 
node. The reactive power into the pth node can be calculated 
from equation (2.11) 



-Im 


Ft* 

' P 




t 

q=1 


■pi <ii 


( 2 . 12 ) 


The voltage must satisfy the relation 

bV ' = ?E , ' 

' P ^ sh , p i 


(2.13) 


in order to calculate 1he reactive node power required to 
provide the sdieduled node voltage. The present estimate 
V* must be adjusted, therefore to satisfy Eq.(2.13). Ihe 
phase angle of the estimated node voltage is 



Im(v|) 

Re(V^) 


(2.14) 


Assuming that the an^es of the estimated and scheduled 

Ic 

voltages are equal, the adjusted estimate for is 


7 


■k 

p(new) 


= Kh,pi ^p + 3Sin 6p) 


(2.15) 


k k 

Substituting in eq. (2.12) the reactive power is 

k+l 

obtained. The new voltage estimate 7^ is calculated 
using Gauss-Seidel iterative scheme, by the following relation 

■yfc+1 Y 7^"^^ - ^ Y 7^ I 

» “ qtl W <1 —1 pq qf 

■ t,. P n ew ^ _j 

(2.16) 
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Similarly new estimates for all the nodes (i=1 , 2, . . . ,n; 

i^S , where S is the slack node) are obtained. 


Ihe rate of convergence can he increased hy applying 
an acceleration factor to the approximate solution obtained 
for eadi iteration. Let a and p be the acceleration factors 
respectively for the real and imaginaiy components of voltage. 
The accelerated value is 


yk + 1 

p (accelerated) 


+ o:Re + JP In ] (2.17) 


and replace the calculated ty V^|2ooelerated) ’ 

entire process is iterated till 


max 

P 


i 





e p=1,2,...,n; p;^S 


(2.18) 


is satisfied. The angles of these converged voltages are 
the desired stable equilibrium angles for postfault state. 


It is to be noted that since all the nodes of the 
reduced system are voltage controlled nodes, the deviations 
in the node reactive powers are insignificant during iterative 
process. Hence there is no need to keep a continuous watch 
on the node reactive po\«/ers to check whether they cross the 
reactive power limit. 

2.5 INSTABLE EQUILIBRIUM POIHTS AHD REGia'f OE STABILITY 
An unstable equilibrium point of Eq.(2.8) may be 

U) 

determined by minimizing by choosing a proper starting 
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point. It can "be readily seen that the determination of 
(2 -1) unstable equilibrium points beccmes increasingly 

unmanageable as systar size (n) increases. In fact this is 
one of the biggest obstacles in the use of Lyapunov stability 
criterion for power systems. However, an accurate determina- 
tion of the Unstable equilibrium point is not necessary from 
the engineering point of view as the value of is not 

very sensitive to small changes in x around the equilibrium 
points [23]. A method proposed by O.L. Varayana [22] known 
as the corner point method takes advantage of this fact in 
the calculation of the region of asymptotic sta!>ility. 


Consider a n-machine power system, I'ig.2.1 shows 
the power system nonlinearity. Por the system considered, 
there are ^^2 ~ say m, nonlinearities in the 

systan. Each nonlinearity represents power transfer 
between two machines, and is a function of one of the rotor 
angular differences. Out of m an@ilar differences only (n-1 ) 
can be selected independently. Therefor e ^ defined by (2.7) 


is partitioned as 



(2.19) 


Each i^L» violate the sector conditions for 

1 , > . *> k. . Let 5x be considered in a (n-1 ) -dimensional 
1 1 1 1 

space such that its axis represents CJ.. Two 



hyperplanes cr^ = ancL<y^ = restrict the space ^ 

around the origin in which satisfies the 

sector conditions. Similarly 2 (n- 1 ) hyperplanes can he 
drawn' to fom a polyhedron, polyhedron, so foijiied, has 

on-1 . j . 

corner points and it is claimed that the minimum value 
of V(x) , obtained by a search among 2^ ^ corner points, yields 
an approximate region of stability [ 23 ], The values of k^: 
and 1^ are determined by th e following relations 

k. 2(6° -6 °) 

i P 

1^ = -Tt- 2(6 p -6 °) : p < q (2.20) 

where (6 p -6 °) is the angular difference of pth and qth 
machines at the stable equilibrium point and it is associated 
wi th i 6 I, 

Por a chosen^ j, a polyhedron can he constructed. 

Since there are ’^0 , say p, ways to choose^p, p poly- 

T] -1 

hedrons can be formed and a search among p.2 corner points 
has to be carried out. Su h an enumeration and computation 
for practical ^sterns is difficult to implement. 

To eliminate the various choices of a machine, 
say r, is selected as reference machine, the choice of linearly 
independent variables of such a way that 

cr. , i =1 1 , 2, . . . , (n-1 ) , involves an angular difference which 
contains the te^mi This results in a set dp, which defines 
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■fche boundaiy of the polyhedron, and has corner points. 

The choice of Ihe corner point which results in a minimum 
value of Lyapunov function among the 2^~ such points is made 
on the basis that each component coiresponds to 
min Ihus, for a machine taken as reference, 

a unique corner point is deteimined. The entire process is 
iterated for all machines taken as reference machine. Hence 
the search is reduced to n corner points. The minimum of the 
V-function evaluated at n such points, gives the region of 
stability. 

The main steps to calculate the region of stabilily 
can be summarised as follows: 



Calculate all the possible angular differences 
( 6^ i = D = 1, 2, ..., n, Md* 


2. Choose a reference machine, say r. 

3. Among all the angular differences of step 1, picte all the 

(n-1) differences (6? -6?), such that i or j = r, hence 

J 

form 0 j. 

4. ' Construct the polyhedron using (2,20). 

5 . Choose a corner point which corresponds to 


6 , 


7. 


min| |k^? , ^ “ ''»2,i.., (n-1). 

Evaluate V(x) if corner point selected in step 5 does not 
violate sector conditions, 

Repeat steps 2 to 6 for all machines taken as reference 
machine. 
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8. Minimum of all values of say 0, obtained in step 6, 

gives the region of stability : 7(x)< 0. 

2.6 LYAPUTOY PUHGTIOW 

n this analysis, a Pyapunov function based on the 
total system energy concept has been used. This function 
was proposed by El-AiDiad and Fagappan [s], but was later 
modified by Willems and Willems [15], because the transfer 
conductance term in it caused the function to have an 
indefinite time derivative around the origin. The function 
is given by the expression 



-0os(6^-6p-(6^-6^-6°+6p Sin (6?-^) 1 

( 2 . 21 ) 

This completes the procedure detailing the application 
of Lyapunov’s direct method to power systems. 

2.7 DYFAMIC EQUIVAIEFOIFG PROOEIXJRE 

'-Lhe method proposed oj Stanton, Podmore and 
Germond[29-3l] is adopted. The development of their equivalent 
is systematic and general, accommodating different types of 
machines and their control equipnents. The basic steps in 
the development are: 

(i) Identification o f coherent generators. 

(ii) Reduction of coherent generator terminal buses. 
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(iii) Dynamic aggregation of coherent generators. 

(iv ) Reduction of load buses. 

2.8 IDMTIPIOATION OP OOHPREITT GPi'fERATORS 

'^his is a key step in the tediniques using coherency 
based dynamic equivalents. 

Definition pf _Ooherent Generators: 

Two generators ’i’ and ^ 3 ' are said to be coherent if 
their angular difference during the transient period" is 
constant to a specified tolerance 

: 6 .(t) - 6 .(t). < e ( 2 . 22 ) 

* cl 

where e is the tolerance. A group of generators is said to 
be coherent if every pair of generators in the group is 
coherent. 

In this thesis a method proposed by Pai and Adgaonkarf333 
is used to identify the coherent generators. This method is an 
improvement over Podmore’ s linear simulation [34 ]. The tedmique 
of identifying the coherent generators is detailed as follows: 

The power system is divided into two parts namely 

(i) study system, in which the system behaviour is of direct 
interest and where a set of fault locations are considered and 

(ii) external system, which is equivalised. The linearised 
swing equations are 
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' d t ' ‘ 


1 dt 


'mi 



d 6^ 
dt 


- ^ca 


i = 1,2, ... ,n 


(2.23) 


In order to eliminate^Pg from equation (2.23), the network 
equations are modified as 



I 







“1 r 

i 5 




n] 






(2.24) 


^ ^ 

where Hqq., Hq^ '^LG ^LL partial derivatives of 

active power with respect to the "bus voltage angles. 

6 = angles at the generator internal hus 
0 = angles at the load buses. 


where 


^^n"b^^n+1,n ^^n+2,n’ * * • ’ ^*^n+k,ni 


^n’ is the numb ec of generators and ’k» the number of buses, 
in equation (2.24), all the angles are measured with respect 
to a generator ’n’ which we consider as the reference generator. 
Prom equation (2Jt4), we have 






h: 


-1 

•IL 


^S-i 


(2.25) 
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by elimlnatir)g dsePg from (2.23) and (2.25) the system equations 
(2,23) are cast in the usual state space form as 

^ ^ S with x(0) = 0 (2.26) 

The introdu ction of a network: fault results in the 
instantaneous application of an accelerating torque to the 
synchronous generators [31 ,34]. The electrical power output 
of the units may be computed by solving the faulted network 
equations with the voltage ^ehind the generator transient 
reactances fixed at the prefault values. Therefore, the 
system models for the faulted and postfault period become 


i = A 2 + 1 J ^(0) = 2 ? 0 ^ ^ (2.27) 

X =. A X ; ~e ’ ^ "^e ( 2 . 28 ) 


where x_ is the state of the system at t = t. and is computed 
from (2.27). Assuming that the system matrix has all the 
distinct eigenvalues, the response of (2,23) is written as a 
linear combination of modes as 


n-2 


X = 


a. X. e 
3=1 ^ ^ 


X.t 


(2.29) 


where 

t; 


= “ . + 3 = eigenvalue of A 

3 3 


X. 

-3 


X*. + jX” = eigenvector corresponding to 
3 3 


a^ = 


a’+Ja^’ 


scalar which depends upon the reciprocal 


basis vectors to X. and initial conditions 
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let there be pairs of complex and (2-n-2— 2m) real eigenvalues* 
The response (2,29) is written as 


M a-jt 

X = H: e Mi; "^ios p.t + b; Si>i p .t) + 


2n-^ 


b .e 


a^t 


j=1 ^ ^ ' 3 ' ■ ^^2m+1 


where 


(2.30) 


t tt __it 

lb • 3. jt 3, * X* 

-D 3 -3 3 "-“3 


- (iP O 


b" = -aV'Xt -a'.Z ". 
-3 3 -3 3 — 3 


1 ^ 3 


/ Y” tt \ 
i > / 


"3 ^3 


Yt_+ 3^^= reciprocal basis vector to 
(X[ + jSp 


In order ter ascertain coherency between the generators 
'p’ and ’q' we need to examine the maximum angular excursion 
between the state variables From (2.30) the 

expression for-6 ^ written as 


^6 


V% = 


® ■ CC -4 1 , ~ ~ 3 

JZ- e ^ f (b* .-b’,)0osP ^t+(b" .-b!' . )Sin P.t- 
L P3 <13 3 P3 93 3J 


2n— 2 C£.t 

(b / - b .) e 3 

a=2m;i 


(2,31) 


where b* • and b" . are the elements of b!, and bV respectively, 
P3 P3 ■”3 —3 

corresponding to the state variable A 6p^* 

Uow the coherency index h^^^ which is a normalised 
bound on the de^Hation of A 6^ in the transient period is 


defined as 
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PI 


-m "' ' 2n=2 ^ '‘ ‘ C2.3^) 

Max [ (lD« ) + tl) r=1,2, ...,n-1 

r j=1 ^3 ro 3 = 2 m +1 ^3 


The generators ’p’ and ' q* are classified as oohererjt if 
^pq < e 


■where e is Ih e toleran ce. is generally between _+2,5. 


2.9 REPUariON- OP OOHERMT GEi^ERATOR BUSES 


The termijoal b'uses of the generators in a coherent 
group are replaced by an equivalent bus ’ t’ . The choice of 

’V * the equivalent bus vol •tag e is flexible, generally it 

"t 

is the average of the individual bus voltages in the group 


i .e . 




1 

n 


n 

E, 

k=m+1 



e. 


1 

T1 


i: 

k=m+1 


9 ,, 


(2.34) 


where 0^ is the angle of voltage 

Eor the purpose of clarity, the pro -^ss of combining 
the terminal buses of the generators in a group y^n be 
interpreted as a series of operations on the physical 
network as explained below. 

(i) Once the voltage Y_^ is defined as in (2.34), 
each generator terminal bus in the group is connected to 
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EEIiAIEDER OP ORIGHAI, NETWOEK 






1 


0 


BOUIEARY BUSES 




^23 


^ I 
& 


5 i " "" 

r 


2-, U 


EBSKIUAIi BUSES 


^0 


0 


PIG .2 *2. a OOUPIGUEAIIOU OP THE COHEEEHT GEIIEEATORS 
IH THE ORIGIHAL SYSTEM 



REMAIHBER OP ORIGnTAL HETWORK 


L 


5 

iar 


6 


o 

5> 


-r * 


4^ 


a 


1t 



PHASE SHIPTIHG 
TEAHSPOEMEES 


PIG . 2. 2 .b THE CCHEEEHT BUSES ARE COIHECTED TO THE EQUIVAIEM} 
BUS THROUGH lEEAB PHASE SHIPTIHG TEAHSPOEMERS. 
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PIG. 2. 


i 


PIG. 2 


REMINDER OP ORIGINAL NEIWOEK 




t 


3.C THE BRANCH BETWEEN THE COHERENT BUSES (2) AND (3) 
IS REPLACED BY THE EQUIVALENT SHUNT ADMITTANCES 
AT BUSES (2) AND (3) . 




REIilAINDER OP 


. IT , V f 



i 

ORIGINAL NETWORK ! 


\ 



. 3. « GENERATION, LOADS AND SHUNT ADETTANCES ARE 

TRAHSPERRED TO BUS t AJJD THE ORIGINAL TERMINAL 
BUSES ARE ELIMINATED, 
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bus ’ t* through an ideal phase shifting transformer with 
complex turns ratio a-j^ =(V^/Y_^), where b is a boundary bus. 
Under coherent conditions, the ratio a^^ is constant and no 
circulating airrent flows through the phase shifting 
transformers. Thus inclusion of the phase shifting 
transformer does not affect the systan (I'ig. 2.2b). 

(ii) In the second step, th e non-radial lines between 
the coherent generators are detected and are to be replaced. 
Since there is no change in power flow over these lines 
during transient period, the lines are replaced by two 
shunt admittances at the respective end buses (i’ig.2.2c) given 


by 


Jbt "" . 



(2.35) 


(iii) The generation, loads, shunt admittances on the 
terminal buses are transformed to the equivalent bus ’t’. 
The shunt admittances are adjusted to account for the 
complex turns ratio of the phase shifting transformers. 

Then the entity of the terminal buses i s removed by 
connecting the phase shifting transformers in series with 
the respective interconnecting lines (Pig. 2, 2d). 

2.10 BITTMIO AGGREGATION OE GMERATING UNITS 

The objective of dynamic aggregation is to determine 
the equivalent generating unit model from the individual 
generating unit models in the group. The parameters of each 
of -Sie components of the equivalent generating unit are 
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obtained by considering the corresponding components of the 
individual units in the coherent group. 


The tedinilue adopted for determining the parameters 
is to numerically adjust the parameters of the ocmponent in 
the equivalent generating unit, so as to obtain a minimal 
error between liie frequency response of the equivalent transfer 
function of the component and sum of frequency responses of 
the invividual transfer functions. 

Since, -the classical model of the synchronous machine 
has been considered the procedure of dynamic aggregation 
will be described only for the rotor dynamics, and synchro- 
nous machine. 


1 , I^pi9r_I>ynamics 


M 


The swing equation 

, ^ _ p — p _ 

i dt ” mi "^ei 


for one machine is given by 

Di^i 


(2.56) 


Since the machines in a coherent group have nearly equal 
speed deviations, the summation of the equations (2.36) over 
the generators in the coherent group yields 

Prom th e equation (2.37), it follows_ that (i) the equi'^alent 
inertia constant is the sum of the individual inertia constants 
in the group and (ii) the equiv^ent damping ooeffi-ient ls 
the sum of the individual damping coefficients. 
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2. Syn iiroTious Machine 

The synohronous machine in a coherent group is represented 
hy the classical model of a constant voltage behind transient 
reactance, hence the terminal voltage is same for each machine 
in the group, since they are connected in parallel to the same 
bus after the reduction of the coherent bus. The equivalent 
reactance is given by the value of the individual reactances 
taken in parallel, 

2.11 REHJITIOW OF LOAD BUSES 

The loads are treated here as constant im; 5 edan ces 
(linear characteristics) hence the load buses can be eliminated 
by G-aussian elimination formula, vhich is also known as 
Ward-Hale or Kronas reduction formula. 


Elimination Formula 

The network equations are written as 


U]] 
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(2.38) 


where i!^ is the total number of buses in the system 
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Operating the matrix with standard Gaussian elimination 
technique for eliminating the variables by columns, at 
i~1 th stage we leach the following structure 


■ 


■11 


'12 


'IF 


'1 ! 


? 0 


0 






I 


(2.39) 
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‘•Fi 




Y 


^^4 


Extraotlag i, i+1 , F rows from (2.39), we obtain 




! ' 


^ y 


ii 


'iF ■ 


Y. 

1 


(2.40) 


T± 

! % 


! 

L 


; 


S Y r 5 

j 


Equation (2.40) is an admittance equivalent of the original 
network as viewed from the nodes i, i+1, F, 

This completes the description of the dynamic 
equivalencing procedure. 


2.12 EESTJLTS OE STUDY 

The system chosen for study as mentioned earlier is 
the UPSEB system which consists of 13 hh chines, 71 buses and 
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a-nd 94 lines. The single liT^e diagram of the sjrstem is given 
in Jig. 2,3. The data of the system is given in Appendix A, 

The results of the dynamic equivalencing -will be 
first prese-nted. The results frcm coherency analysis show 
that generators 9, 11, 12, 13 at buses 44, 64, 66 and 68 
respectively are coherent. 

The bus ITo,44 is now assumed to be the equivalent 
bus and the values of the phase-shifting transformer turns 
fatio and angle are presented in Table 2.1 < 


TABLE 2.1 


ERCM BUS 

-1 mr -i«-W . 

X n _ •% f -- .Ik -1 

TO BUS 

PHASE SHIETING 

'ratio 

TRAHSEORMBR DATA 
ANGLE 

_ ^ » *f . xx; . ^ . 

45 

44 

1.05 

+0.06108 

67 

44 

1.05 

-0.022 

65 

44 

1 .05 

-0.012 

69 

- "X « ^ M 

44 

1 .00 

-0.0296 


The details of the equivalent generator at bus l!ro,44 
are presented in Table 2,2. 

TABLE 2.2 


DETAILS OF 7AIJJE 

EQUIVALENT 

Gil'iERATOR ^ ; ■ 

Total generation 665.0+3130.2451 


Total inertia 11.9530 

To tal reaotan ce 0,0765 
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PIG. 2 . 4 : GENERATORS 9 , 11 , 12,13 iPTSR EQUTVAEENCBrG 
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The change^ in the physical configxiration of the system is 
as shown in Pig. 2.4. 

Six representative fault locations are chosen in the 
study system indicated in Pig. 2,3. The critical clearing 
time for a 3-phase short circuit at these locations is 
calculated hy applying i) lyapuviov* s method to the full 
system, ii) Lyapunov’s method to the reduced system after 
dynamic equivalen ci'^g. The results are compared with the . 
clearing time obtained by the step-by-step base case 
transient stability method. The details are presented in 
Table 2.3. 

TABLE 2.3 


METHOD 


1 ) STEP -BY-STEP 
MSB CASE 
TRAITS lEHT 
STJiBILITI 


2 )LYAPUH0V 
METHOD EOR 
THE PULL 
SYSTM 


3 )LYAPUHO? 
METHOD POR 
THE DYJiAMIG 
EQUI7ALMT 


CLmRETG TM POR PAUL T LO >TIOiTS BT SECS . 


Bus 11 
Line 8 

, . ^ M 1 ■mt. % 

Bus 10 
Line 5 

Bus 9 
Lin e 2 

Bus 17 
Line 2? 

Bus 25 
Lin e 24 

X Ik ‘ xn. ir ‘ ^ x 

Bus 21 
Line 30 

0.23 

0.2 

0.13 

0.16 

0.23 

0.31 

0.13 

0.10 

0.075 

0.095 

0.15 

0.20 

0.13 

0.10 

0.075 

0.095 

0.120 

0.180 








The computation time taken by the DEC-IO system for the three 
methods used in calculating the critical clearing time is given 
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In Table 2,4. The lyapunov method applied to the dynamic 
equivalent results in a considerable saving of computer time. 


TABLE 2,4 



METHOD 

OOMPU TATI Off TIME 

IH SEGS 

i) 

STEP-BY-STEP BASE CASE 

TRAHSIEHT STABILITY 

22 

ii) 

LYAPUiTOV METHOD EOR THE 

TULL SYSTEM 

6,64 

iii) 

LYAPUHOV METHOD EOR THE 

DYHAMI I EQUIYALEi-^f T 

« % jr _ VYV ft * ' — jn r _ \ 

5.45 

2,13 

COHGLUSiaiS 



The Lyapunov function used for the purpose of analysis, 
was constructed by neglecting the transfer conductances. The 
results obtained were slightly conservative compared to the 
actual clearing times 

The dynamic equivalent describes the system quite 
accurately, and is compatible with the Lyapunov’s method. 

The results obtained showed that for fault locations well 
within the study system the clearing times were very accurate 
but for fault locations at the boundary of the study__bystera , 
the clearing times were not the same as those obtained by 
applying Lyapunov’ s method to the fall' system. 



There is a large saving in terms of computation time 
■when the Lyapunov’s method is used and this may he one of the 
facjtors for a trade-off with conservativeness, especially 
when the Lyapunov’s method is used as a preliminaiy screening 
procedure to reduce the number of studies to be carried out 


in detail. 



Siiature is "blie ch^Jiiof a man 
Characier is "the chacnin siaiure 
Enowledge is the chacnin character 
And forgiveness the glory of knowledge, 

•VEDAS’ 

GHAP2ER III 

A KBW DECOMPOSITION STEUCTUEE EOE STABILITY ANALYSIS OE 
LAEOE SCALE POWER SYSTEMS USING VECTOR LYAPUNOV FUNCTIONS. 

3.1 INTRODUCTION 

In this chapter a new decomposition structure for 
stability analysis of large scale power systems using vector 
Lyapanov functions is proposed* Both conceptual and numerical 
difficulties in the analysis of transient stability of multi- 
machine power systems increase rapidly when the size of the 
system becomes large [ 13 , 35]. A new approach to this problem, 
is based on the concept of vector lyapunov ftmctions and the 
decomposition - aggregation technique. This new procedure of 
stability analysis of the complex system is a two level concept 
At the outset the composite system is decomposed into siE^lerr 
subsystems and their interconnections. At the lower level, 
each of the subsystems is tested for its stability* Scalar 
Lyapunov functions are constructed and standard techniques are 
used to estimate stability regions of the low order subsystems » 
On the higher hierarchial level, an aggregate model involving 
a vector Lyapunov function is constructed in order to use 



estimates of subsystem stability regions and determine an 
estimate for the stability region of the overalh system, fhis 
method is suitable for handling large complex power systems 
due to the following properties: 

(i) Ihe method can be carried out systematically to get 
stability region estimates which can be easily interpreted in 
terms of system parameters. 

(ii) In the course of transient stability analysis, decomposition 
can be used to take advantage of (or obtain infoimation about) 
the special structural features of the power systems. 

(iii) The method opens up a real possibility for mo3?e refined 
models of the subsystems to be included in the analysis of 
large power systems. This is based upon the fact that when 
more complete descriptions of the individual machines are used; 
it is easier to construct appropriate Lyapunov functions for 
pairs of machines than to find a single Lyapunov function for 
the entire system. 

Application of the decomposition - aggregation method 
to the transient stability problem of large scale power systems 
was initiated by Pai and Harayana [5t ] . Their results were 
obtained in the right direction tov/ards evolving an effective 
decomposition procedure',' but their technique had one drawback. 
They decomposed a n-machine power system into n(n-1 )/2 subsystems 
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thus making the number of subsystems exceed the order of the 

system whenever the number of machines is larger than three, 

Jocic & Sil^ak [53 land Jocic, Ribbens - Pavclla & Siljak [52 ], 

decomposed a n— machine power system into n—l subsystems, but 

their method suffers from one limitation i,e, a reference 

D h- 

machine n should be so chosen that ^ ^ i = 1 j2, ..,,n-l, 

where D is the damping coefficient and M the inertia constant. 
Thus the decomposition structure hinges on the choice of the 
reference machine. 


In this chapter a nev7 decomposition structure is 
proposed, based on the centre ’ of angle principle originated by 
Stanton [ 54 ] and Tavora and Smith [ 55 1 • 4n appropriate subsystem 
Lyapunov function is constructed using the method outlined in 
[52 ], Using the separate results of ref. [41,45,48,56, 57 I 
stability region estimates are computed by a quadratic Lyapunov 
function. The technique is illustrated for a 4 machine numeri-cal 
example detailed in Ref, [8 


3.2 DECOMPOSIDIOl OR A MULTIMACHINE POvTBR SYSTEM 

The swing equations for the postfault system character^ 
izing a n~machine power system as indicated in equations (2.5) 


of Chapter II are given by 

n 


Ml i'l Ii 5i = - ^ BlE^Ylj [ Cos( S i.-e. . ) _ CCS ( ] 




( 3 , 1 ) 
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Poliowing the method proposed by Stanton [ 541 and 
Tavora and Smith [55] the motion of the generator rotors can 
be resolved into two components: 


(3.2) 


(i) Motion of the * Centre of inertia’ defined by 

g 

where Mm s M. 

^ i=T ^ 

(ii) The relative motion between the centre of inertia and 
individual generators, known as the off-centre angles and 
denoted by 

i = 1 ,2, .«.,n (3.3) 

The off-centre angles are not all independent and have the 
following properties 


5 . « 6 . „ 6 

10 10 




= 0 


n • ' 

2? M. 6.^ = 0 


t Ml'^o = ° 

i io 


(a) 


(b) 


(c) 


(3.4) 


The equations of motion in the new co-ordinates are obtained as 
follows [ 58 ] . 

Prom equation (2.1) of Chapter II 




(3.5) 
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I'rom eq.n, (3.2) 


Mm = fc M. s’. 
° i=:1 ^ ^ 


n . n 

-z:: D. 5, + (p .^ .) 

■ i;=1 ^ i=l 


m. ei" 


Prom the R.H,S, of eqn. (3*6) adding and subtracting 
we get 


6 ^ 5 

o = 5$ ^0 “ 

v/here 


-A- 

i^f 


% 




i=1 





(3.6) 


n 

>' - ^o 
1=1 ^ ° 


(3.7) 


The linear dependence that exists among the off-centre variables 
indicated in (3.4(a)), reveals that it is always possible to 
describe relative motions within the system through n-1 off 
c-entre angles. Manipulation of equations (3.5) and (3.7) yields 
a set of differential equations which explicitly characterizes 
the internal motions of the system. These equations are obtained 
as follows : 


Adding and subtracting 6^ from the E.H.S, of equation (3.5) 
yields 


(3.8) 


Subtracting equation (3.7) from (3.8) yields 

" D. * P. S. P. Pm * 

io “ io Mg, r ^ Mg,^ o 


•fc 


-rr 


^ ^ - P • ^ P . 

f Tn*f- 7n*i- \ / \ n i o ^ ^ 

■ HL ■ ^ -XT - 


(3.9) 



50 


The set of equations (3*7) and (3*9) are equivalont to 
the original system description provided by (3»1) * In 
the new set, however, the equations are expressed in terms 
of variables which have excellent properties for a proper 
casting of the state space form. 


The subsystems are no7\r cast in the state space form. 


To aid in the suitable selection of a subsystem nonlinearity 
"til 

the n machine is taken as the reference machine without loss 
of generality. Using the fact that 6^ = 6^^ “ ^no “ 

ejquation (3.4-(b)), results in equation (3.10(a)), Equation 
(3.10(b)) follows Uirectly from equation (3.9). Because of 
the appearance of in the state space equations, we attach 
equation (3.7) to each of the subsystem descriptions. The 


resulting state space model for each subsystem is given by 

n-1 

-I ..1 V ■' 

5. =6. -6 =0. -cf =ri+M.M 

in io no lo no i n xo n 3 0 ° 


(a) 


D. 




CJ 


rt~1 
N 


'*•- ’H. . V i.,'<w 


>“A ^.1 


JO 


n 




- A - 


„ JL'-o ^ "I “ij *13 

1 X 


n~l n 

1 E C 

i=l ,1=11- 




S' 


[Cos ^ij *■ cos 6 ] 


(b) 
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-1 






n-1 

.-1 


- • - )"lo - W /=T 

3A 

n-1 n 

- 2 S ^id f "i3^ 


where 

=^io. «„=^a.' hd=®i®dhd ®ld-®lj^- 




In (3*10) correspond to the postfault stable equilibrium 

state given by equation (3.1)* Defining an incremental state 
^ .0 

variable the state— vector for each subsystem becomes 


^ ^in" 


T 


in, 


10 , 


(3.11) 


The components of the subsystem state vector are the relative 
rotor angles (any machine being chosen as the reference machine), 
relative velocity between machine i and centre of angle velocity 
and centre of angle velocity itself. This decomposition 
structure differs from that in Refs. [52 53 1 and in particular 

the subsystem has a simpler structure* The interconnection 
terms though different reduce to those of Jocic and Siljak [53 ] 
at the aggregation stage* 

With the choice of the state vector as in (3*11) and 
expanding f^^ = Cos( ^ trigonometrically^ 

equation (3 *10) can be oast in the formi , : 


^£6* N». i*'W 


59478 


imm 
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0 


i 0 ! 




.-1 


Si = 


■ 0 


! 0 


D-- 


D, 


- sc> ! £± •" i -*=1 


% i 


- B 


% 


J 


y. = [ 1 0 0 ] 




i = 1,2,,,,, n-1 


where 

X 


^ ^In" ^ 1 n, " --" ^ 


n-1,n" ^n-1,n, "'-''lO, * * ,o,^o 
is the state vector of dimension (2n--1 ) for the overall system 
( 3 . 12 ) 


(3.12) 


(3.13) 


(y^) = Sin (yj_ + S^) - sin 6 


in 


(3.U) 


forms the subsystem nonlinearity 

i ^•“1^ 

n ..1 




h^(x)=f i -b^lj_cot C ^ ? ^0O 


d=i ^ 


n~l 

^ -u- 9 wr'* • Sos ©. . [Cos 6. .-Cos 5. . ] 


id 


id 


o 

5. . 
3-3 




^Ty.-1 !_/ ^ k. '. Gas ©. . [ Cos 6 . .-Cos 6 . . ] 

— iiivim ^_'| IJ n in Ij 




3-3 


3-3 


3-3 


i -m; 3 »r“n“5“n 5“jo 


-1 


(3.151 


d;^ 
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and 1. = Cos 6 . - Cos 6^. , b. = M7‘ sin © 

j- in ui ' 1 1 m 


m 


iPor a clearer understanding the decomposition is illustrated 
for a 3-machine system, with machine 3 as refei*ence« Ihere will 
be two subsystems given by: 


SUBSYSTEM 1 


i 0 

W 




. . . , >1 ■ ) 




= t 0 


1 


% 


r-1 


i !D-i— B- tM'iM't 

I-o . L ,5„,1-.3. 


y-f - f 1 0 0 ] Si 


0 i 




D-J Dfji f ^ d” ", 


Mj* 


L "J 


(3.16) 


SUBSYSTEM 2 


10 


Xo = ^ 0 -( 


I 

i 

ro 


(1+M2M5^) 

Dg - (I)2-I)^M2M"‘), 
Mo Mm 


0 


t : 

I 1 

i 0 I 




^Ml " j J252(^)+h2(x) 


"T 


D2-B^M2M^'’ 

‘ % ^ " ^ 


Mq>. 


^ 0 i 


12 = [ ■' 0 0 ] X2 


(5.17)! 


This completes the decomposition procedure. 

3.3 SUB SYSTEM AUAIYSIS 

Each free subsystem is given by •eqn.( 3 *18) -obtained. 

from (3,12) by neglecting the interconnection terms. 
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^ 0 




Zi “ n 


(i+Miiqi ) 

0 

"1 

i 

t 

1 

! T 

! 0 i 

> f 


.(h. 



0 % 

i 1 

* 1 


% 

1 

1 

S 

i 

; i 
■ * 

f f 

! °f 

0 0 ] x^ 

• 




(3.18) 


where = (Xj^-j x ^2 ^ 13 ^ state vector defined in (3*9) 

with Xj_^= 6 jjj- 6 ^^ ^13 ="-'o' speoiflea 

in ( 3 * 14 ) • l^he objective of this section is to construct a 
Jjure*— Pos1:nikov type of Lyapunov function using Kalman’s const- 
ruction. The free subsystem (3.18) should be first transformed 
to obtain a subsystem with an asymptotically stable linear part. 

In other words, we add to the right side of the second equation 


1 


in (3.18) the zero term ab^y^l^ -abj_ [0 1 1 where A 
SDTQ ^.^ and a is a positive number, After a simple regrouping 
Qf the terms, we get the transformed subsystem as 


in 


Xf = A^Xj^ + B^>0i(Zi) 

T 

z± = ^i % 


i — 1,2,.,., n— 1 


( 3 . 19 ) 
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where ■ 0 




A -La, 




0 


1 

0 




I » i 


'T i 




0 s- 


(3.20) 


^ J I J 

Kalman’s construction procedure for the above system is detailed 
as follows. 


1 ) The transfer function of the linear part of the system 
under consideration can be written in terms of the parameters 


of 

eqna.(3.19) as 



W( 

s) = (si 


(3.21) 

2) 

let Z(s) = (14-I^s) 

¥(s) 

(3.22) 


so that 




+ Z(*-3‘^)) = 

= Re (1+fj_3^)W(ac^) 

(3.23) 

factorize the expression 

(3.25) as 


1 

S 

(Z(jcj )+Z(- 3 <u!)) = m 

(-jc^)m (a'-O 

(3.24) 

3) 

Determine the scalar 




Y= CO + 


(3.25) 


4) Solve the identity 
m(s) - 

for til© componQii'ts of tiie Tecfor 
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5) Solve the Lyapimov matrix equation 

A? H- + H.A^ s; - 2 -. 2-^ 
a. 3, ^ 1 i ^ 


(3.27) 


for the syMnetric positive definite matrix: H^>0, A point to be 
noted is that is scalar chosen so as to satisfy the Popov's 
frequency criterion 


Re (1+ ) > 0 
The Lyapunov function is given by 

cT X ' 

^l(^i) “ ^i J ^ ^i^3 


j^JiCy^) dy^ 


(3.28) 


where is abtained from in (3.2?) by making ® =0, 
The function the estimates 

*111 


(3^9) 


as shown by Weissenberger [ 56 ] . 


Jocic and Siljak [53 ] have shown that 

ir]^^ = X(H^) ^ j ^2 ~ (Hj^ ) and h i=1 ,2, «. .jn— 1 

(3.30) 

where X and /N are the minimum and maximum eigenvalues of the 
corresponding matrices. is obtained from by setting a=0 
in ( 3 . 27 ) and + i choice of comes 

from the fact that yj.jZ5^(y^) ^ K!_Ji where m is the slope of 
the function /j;(yi) at y^=0, and the fact that the integral in 
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(3*28) can be majorized by ^ ^±y±* 

To derive the second inequality in (3.29) '^i(^i) ( 3 .I 8 ) 
is first computed. A positive number e ^ is so chosen that, 
the inequality 




(3.31) 


is satisfied in an interval 2 c£.j. 1, where the limits 

x^.J and x^.| are the nonzero solutions of the equation 


e sin (x^, + 5^^) - (3.32) 

using (3" *.51)’ is calculated from the inequality 

Vi(Xi)(^^l8) ^ -Xi^ Gi Xj_ (3.32) 

The estimates (3.29) are needed in constructing an 
aggregate model for the system (3.12), The stability region 
estimate for each decoupled system (3.18) can be readily 
computed using the results of Walker and McClamroch [59] and 
Weissenberger [56 ]as 

. V /'■v . i=l,2, .,,,n-1 (3.33) 


where 
V 


X, = V,(xp<-.-V,| 

Here 1 v ^ 

0 . i rA ■ 2 + 

yi=^il»^il a O. O 1 

[Cos 6^ -Cos(yj.t 

and xt^.j I x”j^.j are the solutions of (3.32) • This completes the 
subsystem analysis. 
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3.4 STABILITY EEGION ESTIMATES EOR OVERALL SYSTEM 

An estimate X of the stability region Z for the overall 
system (3.12) is defined. 


A region Z of points Z 6: R*^' where is the m-dimensional 
Euclidean space , is an estimate of the region of 

f 

asymptotic stability Z(ZCZ) for the equilibrium x=0 if and 

only if x=0 is stable in the sense of Lyapunov, (t^jX^) KZ 

where T represents the time inteival, implies xCtit^jX^) ^ Z for 

all t <ST^ where T^ is the semi-infinite time interval (t^, ) 

and lim x(tj t , x ) = 0. 
t ->•<£><» ^ ^ 

To compute Z defined alove a quadjKttic Lyapunov function is used. 


V (x) = 

where the vector Lyapunov function is defined as 

T 

V(x) =[vf (xj), v| (Xg), .... 

and matrix D has the diagonal form 

D = diag d-j > ^2^ *** * *^n— 1 j ’ 
with d^>0 , i=:1,2,..., n— 1 


(3.35) 

(3.36) 


(3.37) 




The desired asymptotic property of Z is established by "tiie __ 
stability of a (n-1 ) x (n-1) matrix ¥ = (w^^-) which is defined 

^ - />"*(%) X(&i) - i = 3 


W, s= 
i3 


(3.38) 




X“<hP A-(Hj)fy 






59 


where ^ is the interconriection eonstanl;, The procediire 
to determine matrices H^, and is given in section 3.5 


The region 


X =fx ; vix) < 

(3.39) 

with ■)/ -s min S d. V. { 

(3.40) 

n-1^ 


and V(x) «= d.Y. (x. ) 

i=1 i i 1 

(3.41) 


is an estimate of the stability region X containing the 
equilibrium xs=Q* of the system described in (3*12), if the matrix 
W defined by equation (3*38) satisfies the ineq'ualities 


1 ^11 

w-j 2^ « • # 

^Ikl 

[ 

^2fc| 

• e ( 

i 

: W21 

> « • • 

W 22 • • • 

K. 

1 ^ki 

^k2 • • • 

1 


> 0 , ^ fcsl • jn^i 


(3*42) 


«• *'s 


To determine the interconnection constants the 

following procedure is adopted. (5,1 2) computed as 

^ , “s . . r Tr 




S ^ [ - (sp lixji ‘ + II 2 fij IFlii (b) 

i='' U± 

(3.43) 


+ [grad ^^(x^)!] h^(x)‘ 
.. 2 


n-1 


(a) 


comparing the two equations in the set (3*43)^ are obtained 


*D© tails givsn in Appendix B 
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Uow using equations (3.34), (3.39) and (3.40) the 
stability region estimate can be computed, This completes the 
procedure for computing an estimate for the stability region 
of the overall multimachine power system. 

3*5 NIMBEICAI EXilPIE 

In this section, an illustration of the proposed approach 
is presented using a realistic four machine power system. The 
data and details of the system are available in ref. |_8 f • 

= 75.35, M 2 =» 1 . 130 , = 2.260, 1^=1.508 


“1 = 

1 , 

^2 = 12 

1)3 = 2.5 D^= 

6.0 

T21-” 

0*66416 

73-1=0'. 66098 

= 0.7606 ■) 
e^1 = 99 . 152 ° ^ 
S41 = 2 . 96 °" j 

Details Eefer to 

®21~ 

, ■» 0 

C2i = 

91.294° 

5 . 85 ° 

0 ^^= 97 . 059 ° 

63^=13.79° 

postfault stable 
equilibrium point 

( 5 . 44 ) 


Machine 1 is chosen as the reference machine. The positive 
numbers t ^ and 'c are chosen as follows : 

=fe^ = 0.5 (3.45) 

Equation (3.32) is solved for y^ = to get the domain for 
subsystem nonlinearities as 

= 99.05°, - - 117.9 

:= 85 . 7 °, 150.074° 

« 103*8°, x^ = r1'l5.34°' 


(3.46 
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which are used to compute V° in (3.34) for the three 
decoupled subsystems. 

The minimum and maximum eigenwalues, which are needed for 
aggregation, are found to be 

= 0.361 , aCh^) = 157.765 , M<^2^ - 1.389 

X(Hj) = 0,00651 , A(l!j) = 567.651 , A(&3) = 0,16517 

X(H^) 0.02399 ,A(H^) = 231.136 , A(G^) = 0.25554 

(3.47) 

By using equation (3,43), the bounds on the no nl inearities 
are calculated as 

1^2 = 0.1627 ?23 = 0.000065 f 24 = 0,000024 

f ^2 = 0,000065 ^33 = 0.0957 f 3^ = 0.00001312 

f^2 = 0.000024 f 43 = 0'.00001 312 = 0. ‘081 9 ( 3.48) 

The aggregate matrix; 

i~0. 00882 0.00134029 0,00025776 I 

W = * 0.00134029 -0.0002906 0.0001005 [ 

1 ■ ■ • ‘ * 

! 0,00025776 0.0001005 -0.001105 | (5.49) 

satisfies the inequalities (3.42) and the matrix D of 
(5.37) with the choice 

dg = dj = = 1 

Prom equation (3.34) 

1° s 29,29, = 0,43758, = 10.16816, X ==0,43758 

(3.5a) 
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Using equations (3.39), (3.40) and (3.41) we get the 
Stability region estimate as 

V(x)= + »^2o( ®21-«21°)- 

+0 . 49 4/-3 . 64i-J, ‘*^20*0. 27 ( Sj-, - 63° ) -’^-3 .64'-’20 --^o 
2 ■ ^ - 
+ 157 . 580 ^ + 0.92 (0.99-(^21'" ^21^°*''^’* " 

)‘'0*057*-^q( ^ 1- ^3i)+i(^3l-^3l)^^0 

+1 .004.u>3/ ~0.057( ^31- ^3?H-17.10 ^30% 

+567.13^§ +0.34 (0.,97 -(651-5^0)0.23 - Cos )+ 

1 .93( ^41 -.^^1 ) 4^ “^4,0*^ '^41“ 4,1 \l ^■^^^4l*"\l ^*^0 

+0,55^4Qr9.67c^^c^40: 0,13(5^1- |°) -9.67<^;o H 

+230.75f-^Q +o;53 (0,99 - C®41-«4^ 0.051 - Cos ^^1) 

< 0.43758 (3.51) 


The faulted system equations are then integrated frcmthe prefault 
equilibrium state, till the region of stability is reached. This 
gives the critical clearing time. The critical clearing time 
obtained for this example is 

t^ = 0,22 sec ( 3 , 52 ) 

as compared to t„ = 0,4244 seconds, obtained by applying 
Lyapunov's direct method as in ref. [8 ]. This completes the 
numerical example for the 4 machine system. Although the results 
are somewhat conBervative, this represents the first application 
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of the vector Lyapunov function approach to a realistic 
power system to compute critical clearing time. Selection of 
parameters in (3.45) was somewhat arbitrary. An optimum choice 
of those parameters will lead to a reduction in the conservativeness* 
This is an open area for research. 

3.6 CONCLUSIONS 

This chapter has demonstrated the applicability of the 
vector Lyapunov functions to the stability analysis of large 
scale power systems. The key to a successful application of the 
method lies in an effective decomposition of the power system. 

In this chapter a new decomposition procedure is developed. 

This enables one to construct Lure’ Postihov type Lyapunov 
functions. Using these functions a scalar Lyapunov function for 
the overall system is obtained via a vector Lyapunov function. 

The critical clearing time is then evaluated. Although the 
clearing time obtained is quite conservative, the possibility of 
reducing conservativeness by including more system details makes 
this new approach an attractive tool for a detailed stability 
analysis of large scale power systems. 



There are no eyes that see like knowledge j 
there are no ;joys as in truthfiilnessi 
there are no pains as in attachment | 
and there is no bliss as in desirelessness. 


‘VEDAS' 


giAPTER IV 

INCIUSIOE OP 10N~KEGIIGIBEE TRAFSPER CORDUCTAITCES IN 
LYAPUNOV PUNCTIONS POR I.!UITIIIA CHINE POWER SYSTEliS 

4.1 INTRODUCTION 

In the transient stability analysis of a power system 
via the direct method of Lyapunov a suitable mathematical model 
of the system is a basic requirement. Por multimachine systems 
this means that, in addition to the governor and exciter 
dynamics, the danger windings in the direct and quadrature 
axis, in the dynamic model of the synchronous machine, be 
taken into accoimt. This necessitates the use of Park’s 
equations, change of reference axis between machine and network 
reference frames etc., resulting in a mathematical model invol- 
ving complicated nonlinearity. Analysis by Lyapunovas method 
therefore becomes practically impossible* Hence, a number of 
simplifying assumptions have hitherto been made so that the 
mathematical model is amenable to analysis by Lyapunov based 
methods. These assumptions have been detailed in section (2.2) 
of Chapter II, 
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Efforts to relax these assumptions have met with 
limited success^ In particular, there have been efforts 
[6,60,14,12] to relax the assimiption regarding the neglect 
of transfer conductances, Ihe transfer conductances are off 
diagonal conductance terms in the admittance matrix reduced 
at the generator intemal nodes. Their magnitude depends on 
the network parameters of the system, location of the fault and 
loading conditions. Although the conductance components take 
small values as compared with sueceptance components,^ it still 
remains an open problem whether or not we can neglect the 
conductance components from the standpoint of the transient 
stability, Eor some cases, the neglect of transfer conductances 
may be justified, as shown by Ribbens-Pavella [12 ], However, 
if the system is heavily loaded, it has been pointed out by 
Uyemura et al [60 ] that there exists a danger of judging a 
practically unstable system to be stable if we use a l 3 ?apunov 
function which neglects the transfer conductances. 

In this chapter a method proposed by Pai and Varwandkar 
[62] has been used to construct Lyapunov functions which include 
effects of transfer conductances. Computations Of 03?itieal 
clearing time for a practical 13 machine 7^ 6us powcy system 
are done'^ and compared with those obtained in Chapter II by 
neglecting the transfer conductances. 
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4,2 miATEEMATIGAL MODEL AUD LYAPTMOY EUWCTIONS 

Ilhe model of a n—machine power system with, transfer 
conductances can be written as 

X = As - B^f( a ) - B 2 g( <? ) (4.1) 

<7 = Gbc 

where A^B^, B2 and C are constant matrices, x is a (2n-1)- 
dimensional state vector and the n(n-1)/2- dimensional non- 
linearity vectors f ( cr) and g( o ) are such that their i^^ 
component depends on only the i component of <^51.6, 
f^( <J) = and g^( o) = g^( 0^). 


Without transfer conduct^ces, however, the model becomes 


X = As - B-|f( <7 ) 
^ = Ox 


(4.2) 


Typically, for a threeHnachine system with machine 1 as 


reference, we have 
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II 

® 1^2^12 

[ Sin 

( + 

®12^ 

-- 

Sin ^ -j 2 

] 

fjC-’j) = 

E.j E^B-j ^ 

[ Sin 

+ 

CM 

b 


- 

Sin 

] 

f3( ''j) = 

^ 2 ^ 3^23 

Csin 

( O' 

«23^ 

- 

Sin 62^ 

1 


0 

-1 


and similarly 


g-j ( O' •] ) = ® 1^2^12 1"^ ^12^ *” 6 ■J2 1 

S 2 (^ 2 ^ = ^^2 +^ 13 ^ - 

gjC <^5) = 112^5*^23 ^ '^3'^ *^23^ " ^23 ^ 


0 


1 ^ 

* 

t 

^ 4 


(4.3) 


where "tlie double scrip "bed B*s and G*s are "tlie "transfer 
susceptance and transfer conductance terms, respectively. 

The nonlinearity f^( o'jj^) satisfies the sector condition 


0 < 


- <r. 


< K , 


K =5 00 


(4.4) 


in the interval 
0 


- 26 2 ^pq 


^ % 


(4.5) 
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where p and q (q> p) are the indices of the angles at the 
generator internal nodes on which the i^^ nonlinearity depends. 


The nonlinearity g^( unfortunately, does not 

satisfy a similar sector condition. One can easily establish 
that near the origin it is in the first and third quadrants 
if Gpq Sin 0 "biie second and fourth quadrants 

if Gpq^ Sin ^ 0 P* Because of this uncertainty 

in the sign of this nonlinearity, the system (4*1) cannot be 
cast in the Iiire* foim, and therefore the systematic procedure 
given by Pai and Murthy [12] and Willems [26 ] cannot be applied: 
to form a Lyapunov function. 


Now consider the nonlinearity h^( ^j_)=f 3 _( ^ 

"p“<irp(i pc^'pq' ' p^’^pq^ 


G, 


PI 


= E E Y [Sin( ° +© )- Sin ( 6° '+© ')] 

n dicq ^ r occ na'^ '* ncD na''-' 

Y„ „ Sin ©.. 


PI: 


pq 


B, 


PI 


= Oos 

pq. p<i^ 


(4.6) 


This nonlinearity satisfies the Popov sector condition in the 
interval. 


- ^--2( 6°^^p^) ^ Oi < -2(5^^-.^pq^) 


(4*7) 


Introducing this nonlinearity, system (4*1) osn Bg rewritten as 

X = - B^h(d) - [Bg - ] £( ) 

cr = Cbc C4n»8) 

This particular representation has certain advantages 
over (4.1) and (4.2) in that the nonlinearity hj_( o^) contains 
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"transfoi* conducbance "fcerms and still satisfies the sector 
condition* Also, it is possible to show that it is the last 
term i*e, g( ) in (4.8) which presents difficulty 

in forming the Lyapunov function, and if neglected, the 
Lyapunov function can easily be constructed. 

The Lyapunov function for (4.2), i.e. the system 

neglecting transfer conductance^' !_12] , [I 5 ] is given by 

Cbc 

V(x, «^ ) = x^Px + 2 f^( cr)Qdcy (4.9) 

o 

where Q = identity matrix, with a negative semidefinite 
derivative along (4.2) given by 

V(x,cr ) =t ~x^ LL^x (4.10) 

where matrices L & P (P ^ 0) satisfy the equations 

A^P + PA ;= -LL^ 

and PB^ = A^ 0® (Q (4.11) 

The validity of the Lyapunov function (4.9) for the model 
(4.1) can be considered by differentiating (4.9) along (4.1). 
Recognizing the fact that CB2=:0 

V(x, a) « ^(x, o) - 2x^PB2g( cr ) (4.12) 

It is possible to show via a three-machine example 
that certain terms in (4.12) cannot be integrated to form a 
suitable Lyapunov function, Por a three-machine ssrstem [I 5 ] 
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0 

0 

I • 

0 

0 


0 0 

M2_ 0 

0 Mg 

• » • 

0 0 

0 0 


0 

0 

0 

I I 

0 

0 


0 

0 

0 

0 

0 


(4.13) 


so that 


V(x, cy)«v(x, cr )-2 (■^i + 4)gi ( )+(^'l+ <^ 3 )g 2 ( 

Since =?<^-| -^ 2 > ^2 ^3 ~^2 ■“^ 3 » 

V(x, or) = V(x, a)- 2g^ ( a^)or^ - 2 g 2 ( < 72)52 - 2 g 3 ( ^ 3 ) <^3 

“ 4 [<^ 2 Si ■[s 2 ( '^ 2 ^'‘S 3 ( ^ 3 )^ 40 ^] ( 4 . 14 ) 

* 

In ( 4 . 14 ) all terms except V(x, a ) are sign indefinite. Hence, 

almost eyery research work in this area so far has preferred to 

ignore the terms containing g( <7 ) in (4.1). However, the second, 

third and fourth terms in (4,l4i.), although sign indefinite, can 

be integrated and subtracted from (4.9) to give 

Ox Ox 

^ + 2 ^ :^( d )^ o +2 / g®( a )^ a 

Ox 


«x^Px + 2^/ h (a)iQior (4.15) 

which is positive definite around the origin in view of (4.6) 
and (4.7). 

The last term in ( 4 . 14 ), which cannot be integrated 
directly is, in fact, the expanded form of 

2x% [Bg-B., 1 g( a) (4-16) 
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Since (4«8) is equivalent to (4*1) and noting the structure 
(4*1 6) j it can be asserted that (4*^5) constitutes a valid 
Lyapunov function for the approximate model obtained from (4«8) 
by dele .xng the term [B 2 — g( cf), i,e, for the system 

X = Ax - B^h( cr ) 

= Gx ( 4 . 17 ) 

which is in the Lure * form, 

Lq ( 4 * 17 ) all the transfer conductances appear ±n 
h( e), and hence the critical clearing time computed on the 
basis of the Lyapunov function (4,15) is expected to be closer 
to that obtained by the actual simulation of the swing curves « 
This': completes the mathematical modelling and construction of 
the Lyapunov function. 

4.3 KESULTS OB STUDY 

The study is carried out on the system detailed in 
Appendix A, The same set of representative fault locations, 
as those used in Chapter II are taken to obtain an effective 
comparison# Three phase faults are considered at the above 
mentioned locations and critical clearing time is calculated 
using Lyapunovas method for the model including transfer 
conductances and the results are compared with those obtained 
in Chapter II by neglecting ±he transfer conductances# Table 
4#1 effectively summarises the results and it is seen that the 
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es'bimatG obtained by* including tTansfer conductances is 
slightly on the higher side and is much closer to the actual 
value obtained by the step-by-step base case transient stability 
method* 

TABLE 4>1 

METHOB Bus Yf Bus 10 Bus ' 9 ' Bus* 1? to** 2*5 'to' ■ -'21 

Line 8 Line 5 Line 2 Line 27 Line 24 Line 30 


1) STEP-BY-STEP 
■BASE CASE 
TEAUSIEBT 
STABILITY 



0.13 0,16 0,23 0*31 


2.)LYAPUH0V 

.METHOD • • • . . . - 

IIEGLECTIHG- 0.130 0,100 0.075 0.095 0.130 0.20 

TBANSPEE 

CONDUCTANCES .......... 


3)LYAPUN0Y 

METHOD ■ ■ ■ ' • ' ■ * ■ * 

INCLUDING 0.145 0.120 0.085 0.105 0,145 0.21 

TEANSPEE 
CONDUCTANCES 



4.4 CONCLUSION 


In this chapter the problem of transfer conductances 
in Lyapunov functions for multimachine power systems has been 
discussed. A method of accounting the transfer conductances, 
proposed by Pai and Varwandkar [62] has been used, on a realistic 
system. The results obtained, show that including the transfer 
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conductances, gives higher clearing times, closer to the 
actual values, thus establishing the superiority of the 
Lyapunov function based on the model ( 4 * 17 )* 



Man makes the cosmos and its construction the pivot 
of his emotional life in order to' find in this way 
the peace and serenity which he cannot find in the 
narrow whirlpool of personal experience... The 
supreme task • • » is to arrive at those universal 
elementary laws from which the cosmos can he huilt 
''^P by pure deduction. There is no logical path to 
those laws 5 only intut^on resting on sympathetic 
understanding of experience can reach them ... 

’EIHSTEIN* 


CHAPTER V 
COHCLRSIOHS 

5.1 SIMMART AND CONCLUSIONS 

The development of Lyapunov theory for the stability 
analysis of power systems has been carried out with sufficient 
interest in recent years. The method has proven to be an 
effective alternative to the simulation techniques both 
practically and theoretically. Considerable interest, has been 
shown by research workers as far as the construction of Lyapunov 
fimctions is concerned j but the application of the method to 
realistic power systems has been slow. Apart from the inherent 
conservativeness of the method'y the size of the system and 
detailed modelling have been the major stumbling blocks in the 
successful application of this method. Hence, alternate 
approaches are called for to bypass arid eliminate these draw- 
backs. Vector Lyapunov functions and dynamic equivalencing are 



"bwo of 'th.e most promising solutions for the above mentioned 
problems. 

In this thesis, stress has been laid on some aspects 
of the above two problems, The following are some of the 
conclusions arrived at. 

In Chapter II, a coherency based dynamic equivalent 
is obtained for a practical electric utility in India, Iyapunov*s 
method is applied to the equivalent and critical clearing time 
is calculated for various fault locations, Ihe results are 
compared with those obtained by applying Lyapunov's method to 
the full system. Although the results may bo conservative, the 
dynamic equivalent gives clearing times which tally with the 
actual system everywhere except at the boundary of the study 
system. It is also to be noted that compared to using Lyapimov’s 
method for the full system the method using dynamic equivalents 
results in a 20 ^ saving in computation time. This saving will 
be much more for bigger systems. Thus the method is attractive 
for planning purposes and transient security analysis using 
Lyapunov’s method. 

Chapter III presents a new decomposition structure, 
based on the centre of angle principle, for vector Lyapunov 
functions, Ihe vector Lyapimov function approach is a two 
level concept and utilizes the decomposition-aggregation technique 
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for the stability analysis of large scale power systems'. 
Decomposition of the system is a key step in the application 
of those concepts, The decomposition is pairwise in nature, 
hence it is possible to include system refinements in each of 
these subsystems at the lower level. Studies were carried out 

on a realistic four machine system and even though the critical 
clearing time obtained was conservative, the possibility of 
including system refinements clearly outweighs this disadvantage. 
Ways of reducing the conservative nature have been pointed. 

The presence of transfer conductances in multimachine 
power systems and how to accommodate it in the Lyapunov function 
dis the topic of Chapter IV , This fits into our philosophy 
of trying to make the Lyapunov method an effective practical 
tool',' by taking into account a more detailed representation of 
the power system. An. approximate model including the effects 
of transfer conductances is used. Calculation of critical fault 
clearing times in a realistic 13 machine 7 I hus system demonstrates 
the advantage of including transfer conductances, 

5,2 SUGGESTIONS FOR TiJTUEE EESEiRCH 

On the basis of the investigation carried out in this 
thesis, some of the imsolved problems that need further 
attention are: 

1, In Chapter II further investigation is necessary regarding 
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inclusion of voliago regulator dynamics j as well as other 
refinements laJte saturation, in the lyapunoy analysis* 

a* The decomposition achieved in Chapter III usGs.cmly the 
classical model. An important feature of the new method is its 
ability to handle complexity of the power systems, This opens 
up a real possibility for introducing more refined mathematical 
models of multimachine systems. Also one can examine newer 
decomposition methods such as those using Diakoptic techniques, 

3* Computer algorithms can be used to determine the parameters 
'f^f ^i provide an optimal choice of the parameters and e 
with respect to the size of the estimate Z for the overall 
stability region. 

4,* Inclusion of transfer conductances, in an exact manner, in 
the mathematical model for obtaining suitable lyapimov functions 
for power systems is still open for further research, ^e model 
used in Chapter IV, should be exploited further for modifying 
the Lure* type lyapunov funetd-on with the deleted terms taken 
as perturbations in the model. 
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APPEBUIX a 

DATA FOR 13 M/C, 71 gus SYSTEM 

Number of busses ;= 7I Number of machines := I3 

Number of lines 94 Number of shunt loads = 23 


Base power (MVA) = 200 


Bus 

No. 

Name 

Generation 

Load Power 

Bus voit^e s 
Maff. Ans 

1 

OBRAEXTIY 

630.6672 

169.7222 

0.0000 0.00001.0300 

0.000 

2 

OBRAEXT4OO 

0.0 

.0.0 

0.0 

0,0 

1.0577 

-5.931 

3 

0BRA(T)11 

506.0000 

149.4929 

0,0 

0,0 

1,0250 

-1 ,406 

4 

0BRA(T)220 

0.0 

0,0 

0,0 

0,0 

1,0545 

-6,061 

5 

0BRA(T)132 

0,0 

0,0 

0,0 

0,0 

1.0476 

-5,C^4 

6 

0BRA(H)11 

98.0000 

31 ,9867 

0,0 

0,0 

1 .0250 

0,397 

7 

0BRA(H)132 

0.0 

0,0 

12.8 

8.3 

1^0458 

r5',723 

8 

RIH11 

297,0 

124,2276 

0,0 

0.0 

1.0250 

-0.204 

9 

RIHI32 

0.0 

o;'o 

185.0 

130.0 

1.0436 

-4.976 

10 

ROBI32 

0.0 

o;o 

8o7o 

5o;o 

1,0268 

-7.413 

» - - 

11 

MUG1 32 

o‘.o 

o.o 

155'.0 

96,0 

1^0237 

-11^762 

12: 

MUG220 

o.o 

0^0 

olo 

o;o 

1.^0030 

-10.057 

13 

MUI32 

0.0 

0^0 

0 

0 

0 

62.0 

1^0072 

-13;‘829 

14 

mU220 

0.0 

0.0 

0,0 

0,0 

1.0041 

-11984 

15 

G0R11 

184.08 

113.0 

0.0 

0.0 

1.0055 

-9.413 

16 

GORI32 

0.0 

0*0 

73.0 

45.5 

i;0248 

-16.682 

17 

G0R220 

0.0 

o;o 

36,0 

22.4 

1,0051 

-14.460 
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f -t- 


Bus 

Bo. Bame 

Generation 

load 

power 

Bus Voltages'.’ 







Mag. 

Ang, 

18 


0.0 

0,0 

16.0 

9,0 

1 .0097 

-17.881 

19 

BAST 132 

0,0 

0.0 

32.0 

19.8 

0.9812 

-19.786 

20 

GOBI 32 

0.0 

0.0 

27,0 

16.8 

1.0030 

-21.858 

21 

BAH 32 

0.0 

0,0 

32,0 

19.8 

1.0168 

-21 .376 

22 

SUB220 

0.0 

0.0 

o',o 

0.^0 

1,0082 

-17^625 

23 

SUIil 32 

0.0 

0.0 

75 '.0 

46 '.6 

0 ^ 

1i0361 

-19.815 

24 

SUI400 

0*0 

* 

oi'o 

o'.o 

o'.o 

• . - 

o;'9859 

-16^193 

25 

AIL1 32 

0.0 

l.W 

0.0 

* 

133.0 

82.5 

1 .0391 

-17.441 

26 

ALL220 

0.0 

0.0 

0.0 

0.0 

1,0160 

-14.382 

27 

MAT1 1 

304.0 

76.2874 

0,0 

0.0 

1.0250 

-5.453 

28 

MAT400 

0,0 

0,0 

30.0 

20,0 

1.0522 

-12.342 

29 

pam:ii 

261.0 

70.5069 

0.0 

0,0 

1 i0250 

-151519 

30 

PABK1 32 

0.0 

0,0 

120.0 

74.5 

1 .0574 

-19’.947 

31 

PABK220 

0.0 

0,0 

160.0 

'99.4 

1.0148 

-22.817 

32 

pabe:400 

0.0 

0,0 

0,0 

0,0 

1 ,0252 

-15.936 

33 

LUGK220 

0,0 

0,0 

0.0 

0,0 

0,9984 

- 25,515 

34, 

iiU0Ki3a 

0,0 

0,0 

112,0 

69.5 

1,3066 

-29.639 

35 

I.UCK400 

0,0 

0,0 

0.0 

0,0 

. o‘.9664 

-21.803 

36 

SIT 132. 

0,0 

0,0 

5o;o 

32.0 

0.9943 

- 34.050 

37 

SHA132 

0,0 

o';o 

147‘.0 

92.0 

o;9673 

-37.982 

38 

BARI 32 

0,0 

0^0 

93^5 

88.0 

0^9725 

-37.638 

39 

KHAT11 

25".0 

30t3833 

V ,, 

0^0 

o'.o 

1 .0250 

—34.166 

40 

KHAT1 32 

0,0 

o;o 

0^0 

0.0 

1'.'0183 

-36:913 
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Bus 

Fo, Name 

G-eneration 

Load 

power 

Bus voltages ' 
Mag, Ang. 

41 

M0R1 32 

0,0 

0,0 

255.0 

123.0 

1 .0074 -34f.778 

42 

M0R220 

0*0 

0.0 

0.0 

0,0 

0,9791 -31,481 

43 

M0E400 

0*0 

o*'o 

0,0 

0.0 

0,9462 -27.835 

44 

EAM11 

180.0 

55.0370 

0,0 

0.0,0 

1,0250 4.24;448 

45 

EMI 32 

0,0 

0,0 

0.0 

0.0 

1.0483 r30;5l4 

46 

KE1H132 

0,0 

o'.o 

78.0 

38.6 

1 .0200 -32.678 

47 

MI220 

0.0 

o';o 

234.0 

145.0 

o'. 9893 -29.053 

48 

HAEJ11 

341 '.0 

256.0 

o'.o 

0,"0 

1.0051 -25‘;78l 

49 

HAEJ220 

o'^o 

o'.o 

• 

295.0 

183'.‘0 

- ♦ 

o';9972 -30.'500 

50 

MUE220 

0.0 

o."o 

40'.0 

24 ':6 
• 

0.9776 -32';056 

51 

MUE1 32 

0.0 

0.0 

227.0 

142,0 

1,0051 -354309 

52 

MUE400 

0,0 

0.0 

0.0 

0,0 

0.9579 -26.552 

53 

I/IEEE220 

0.0 

0.0 

0.0 

0.0 

0.9723 -32,845 

54 

1EER132 

0.0 

0,0 

108.0 

68.0 

1.0049 -35.318 

55 ! 

SHAM220 

0,0 

0,0 

25.5 

48.0 

0.9852 -31 ".31 3 

56 

SAHA220 

0.0 

0,0 

0,0 

0.0 

r,0l30 -29.253 

57 

SAHA1 32 

0.0 

o';o 

55 ‘.6 

35.6 

1,0162 -29,967 

5J8 

HAEDlSa 

0.0 

0.0 

42.0 

27.0 

1,Ol84 -30,261 

59 

E00132 

0.0 

0.0 

57^0 

27,4 

1.0135 -30,715 

60 

E00220 

0,0 

0,0 

0.0 

0.0 

1,0087 -29.974 

61 

EIS220 

0,0 

0.0 

0^0 

0,0 

1,0197 -274874 

6a 

EISl 32 

0.0 

0,0 

40.0 

27,0 

1.0443 -28.940 

63 

DEH132 

0.0 

0.0 

33^.2 

20,6 

1'.0426 -27'i695 



37 


Bus 

Wo. 

Wame 


GexiGration Load power 

Bus voltages 
Mag. Ang. 

64 

YAMII1 1 


300,0 

72.8823 0.0 

0,0 

1,0250 -19.401 

65 

YAMII220 

0.0 

0.0 0.0 

0,0 

1.0572 -24.68 

66 

lAMIjIVII 

96.0 

25.6366 0,0 

0.0 

1 ,0250 -20,610 

67 

YAMI,I713^ 

0,0 

0,0 14.0 

6.5 

1.0556 -25.641 

68 

MAW11 


89.0 

26.6892 0.0 

0,0 

1,0250 -19.829 

69 

MAW220 


0.0 

0.0 0,0 

0.0 

1 .0508 -25,563 

70 

HABBI 32 

0.0 

0.0 l-t.4 

7.0 

0,9987 -34,790 

71 

SEA220 


0.0 

0.0 0.0 

0^0 

0,9274 •r35,329 





GEWERAPOR BATA 


Gen .Wo. Bus Brom Wame 

Inertia Trans, React, 

w ^ III _ 1 * t. m- V V r 1 - » 

Bamp Const’, 

1 


1 

0BRAEXT1 1 

31 .2500 

0.0432 

O.'O 

- * * 

2 


3 

0BRA(T)11 

15.4050 

0.0567 

O’.O 

3: 


6 

0BRA(H)11 

1.9230 

0.5278 

0.0 

4 


8 

RIH11 

6.6480 

0,2010 

0.0 

5 


15 

G0R11 

6,2500 

0,2160 

0.0 

6 


27 

MAT1 1 

12.3750 

0.1310 

0.0 

7 


29 

PAWX11 

9,2500 

0.1525 

0,0 

8 


39 

KHAT11 

1 .0350 

1 .0821 

0.0 

9 


44 

RAMlI 

2.2590 

0,3520 

0*0 

10 


48 

HARJ11 

12.6710 

0,0648 

0,0 

11 


64 

YAMIIll 

6.1880 

0,15jOO 

0^0 

12 


66 

YAMIjIVll 

2,2380 

0.5280 

0.0 

13,' 


68 

EA.W11 

1,2480 

0,6000 

0,0 
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IiIME DAIA 


iiine 

Ho . 

'‘'TronT 

Bus 

Bus 

line Impedance 

' " Xy ‘ ' " 

^ charge 

Turns ! 

^ A*. .AT 




*■ .A, 

.. . 


1 

9 

8 

0.0 

i * ^ 

0.0570 

0.0 

1 .^05 


9 

7 

0,0520 

0.0780 

0.0090 

1.00 

3 

9 

5 

0.0660 

0.1600 

t 

0.0047 

1.00 

4 

9 

10 

0.0520 

0.1270 

0,0140 

1.00 

5 

10 

11 

0.0660 

0.1610 

0.0180 

1,00 

6 

7 

10 

0,0270 

0.0700 

0.0070 

1.00 

7 

12 

11 

0.0 

0,0530 

o‘,o 

0.95 

8 

11 

13 

0.0600 

0,1480 

0.0300 

1 ,00 

9 

14 

13 

0.0 

0,0800 

0.0 

1 ,00 

10 

13 

16 

0,0970 

0.2380 

o', 0270 

1.00 

11 

17 

15 

0,0 

0,0920 

0,0 

1,05 

12 

7 

6 

0.0 

0.2220 

0,0 

1 ,05 

13 

7 

4 

0,0 

0,0800 

0.0 

1.00 

14 

4 

5 

0.0 

0.0330 

0>0 

1,05 

15 

4 . 

5 

0.0 

0,'1600 

0,0 

1.00 

a 

16 

4 

12 

0.01 60 

0.0790 

o', 0710 

i':oo 

1? 

12 

14 

o';oi6o 

# 

0'.'0790 

• r / 

o'. 0710 

I'^'oo 

* V * 

18 

17 

16 

o'.'o 

0.0800 

o'.‘o 

0.95 

19 

2 

4 

0.0 

0.0620 

0.0 

1.00 

20 

4 

26 

0.0190 

0.0950 

0.1930 

1.00 

21 

2 

1 

0,0 

0,0340 

0.0 

1.05 

22 

31 

26 

0,0340 

0.1670 

0,1500 

I.OO 



Lino 

Wo. 

“'“'Prom 

Bus 

“ Lo “ 
Bus 

line Impedance 

® charge 

Turns ratio 

23 

26 

25 

0.0 • 

0.0800 

0,0 

0.95 

24 

25 

23 

0.2040 

0,5200 

0,0130 

1.00 

25 

22 

23 

0.0 

0.0800 

0,0 

0.95 

26 

24 

22 

0.0 

0,0840 

0.0 

0.95 

27 

22 

17 

0,0480 

0.2500 

0.0505 

1 .00 

28 

2 

24 

0.0100 

0,1020 

0.3353 

1 .00 

29 

23 

21 

0,0366 

0,1412 

0.0140 

1.00 

30 

21 

20 

0,0720 

0,1860 

0.0050 

1.00 

31 

20 

19 

0.1460 

0.3740 

O.OlOO 

1,00 

32 

19’ 

18 

0.0590 

0.1500 

0.0040 

1.00 

33 

18 

16 

0.0300 

0.0755 

0,00801 

1.00 

34 

28 

27 

0.0 

0,0810 

0.0 

1.05 

35 

30 

29 

0.0 

0,0610 

0.0 

1 .05 

36 

32 

31 

0.0 

0,0930 

0.0 

o'. 95 

37 

31 

30 

0.0 

0.0800 

0.0 

0,95 

38 

28 

32 

0.0051 

0,0510 

0.6706 

1.00 

39 

31 

33 

0.0130 

0,0640 

0.0580 

1 .00 

40 

31 

47 

0.0110 

0.0790 

0.1770 

1.00 

41 

2 

32 

0,0158 

0.1570 

0.5100 

1,00 

42 

33 

34 

0.0 

0,0800 

0.0 

0^95 

43 

35 

33 

0,0 

0,0840 

0.0 

0.95 

44 

35 

24 

0.0062 

0.0612 

0.2012 

1,00 

45 

34 

36 

0t0790 

0,2010 

0.0220 

i.'oo 
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iiine 

Fo. 

Prom 

Bus 

To 

Bus 

Lins ; 

Empodanoo 

'^"^chargo 

Turns ratio 

46 

36 

37 

0.1690 

0,4310 

0,0110 

1.0Q 

47 

37 

38 

0.0840 

0.1880 

0,0210 

1 ,00 

48 

4,0 

39 

0,0 

0,3800 

0,0 

1 ,05 

49 

4tO 

38 

0,0890 

0.2170 

0,0250 

1.00 

5.0 

38 

41 

0,1090 

0,1960 

0,0220 

1.00 

51 

41 

51 

0.2350 

0.6000 

0.0160 

1 .00 

52 

42 

41 

0,0 

0.0530 

0,0 

0,95 

53 

45 

42 

0.0 

0.0840 

0.0 

« 

0.95 

54 

47 

49 

0.0210 

0,1030 

0.0920 

1 .00 

55 

49 

48 

0.0 

0.0460 

« 

0,0 

1 .05 

56 

49 

50 

0,0170 

0.0840 

* 

0.0760 

• 

1.00 

57 

49 

42 

0,0370 

0.1950 

0.0390 

1,00 

58 

50 

51 

0.0 

0.0530 

0.0 

0.95 

59 

52 

50 

0,0 

0.0840 

0.0 

0.95 

60 

50 

55 

0.0290 

0,1520 

0,0300 

1,00 

61 

50 

53 

0,0100 

0,0520 

0,0390 

1 ,00 

62 

53 

54 

0.0 

0,0800 

0.0 

0.95 

63 

51 

54 

0,0220 

0.0540 

0,0060 

1,00 

64 

55 

56 

0,0160 

0.<B50 

0.0170 

1.00 

65 

56 

57 

0.0 

0.0800 

0.0 

1 .00 

66 

57 

59 

0.0280 

0.0720 

0,0070 

1.00 

67 

59 

58 

0,0480 

0,1240 

0.0120 

1,00 

68 

60 

59 

0.0 

o.ceoo 

0,0 

1,00 
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Line 

No. 

From 

Bus 

To 

Bus 

Line Impedance 

® charge 

Turns ratiou 

69 

53 

60 

0.0360 0,1840 

0.0370 

1 .00 

70 

45 

44 

0.0 0,1200 

0.0 

1,05 

71 

45 

46 

0,0370 0.0900 

0,0100 

I* 00 

72 

46 

41 

0,0830 0,1540 

0,0170 

1:,00 

73 

46 

59 

0,1070 0.1970 

0,0210 

0 

0 

74 

60 

61 

0,0160 0,0830 

0,0160 

0 

0 

75 

61 

62 

0.0 0,0800 

0.0 

0':95 

76 

58 

62 

0,0420 0.1080 

0.0020 

•m : 

0 

0 

77 

62 

63 

0.0350 0.0890 

0.0090 

1 .00 

78 

69 

68 

0.0 0.2220 

0.0 

1'.05 

79 

69 

61 

0.0230 0,1160 

0.1040 

1 .00 

80 

67 

66 

0.0 0,1880 

0.0 

1 .05 

81 

65 

64 

■0,0 0.0630 
* . « 

0.0 

1,05 

8a 

65 

56 

0,0280 0.1440 

0,0290 

1 .00 

83 

65 

61 

0.0230 0,1140 

0.0240 

1.00 

84 

65 

67 

0,0240 0,0600 

0.0950 

1 .00 

85 

67 

63 

0.0390 0,0990 

0.0100 

1 .00 

86 

61 

42 

0.0230 0,2293 

0,0695 

1 ,00 

87 

57 

67 

0,0550 0,2910 

0,0070 

1.00 

88 

45 

70 

0,1840 0,4680 

0,0120 

1.00 

89 

70 

38 

0.1650 0,4220 

0,0110 

1 .00 

90 

33 

71 

0, 05^70 0,2960 

0,0590 

1.00 

91 

71 

37 

0.0 0,0800 

0 

' 

0 

0.95 

92 

45 

41 

0,1530 0,3880 

0,0100 

1,90 

93 

35 

43 

o', 01 31 0,1306 

0,4293 

1,00 

94 

32 

52 

0^0164 0^1632 

0.5360 

1.00 
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SHUWT LOAD DATA 


S.No. 

Bus No. 

Shunt load 

Admittance 

1 

2. 

0.0 

■ - 0.4275 

2 

13 

0.0 

0.1500 

3 

20 

0,0 

0,0800 

4 - 

24 

0.0 

-0:,2700 

5 

28 

Q.O 

rO *3375 

6 

31 

0,0 

0.200 

7 

32 

0.0 

-0,8700 

8 

34 

0,0 

0.2250 

9 

35 

0.0 

-0.3220 

♦ 

10 

36 

0.0 

* 

0.1000 

11 

37 ■ 

0.0 

¥ 

0.3500 

12 

38 

O 

O 

0.2000 

13 

41 

0,0 

OiOOOO 

U 

43 

0.0 

-0','2170 

15 

46 

0.0 

0.1000 

16 

47 

0.0 

0.3000 

1? 

5 J 0 

0.0 

0,1000 

18 

51 

0,0 

0,1750 

19 

52 

0,0 

-0^2700 

20 

54 

0.0 

0,1500 

21 

57 

0.0 

0,1000 

22 

59 

0.0 

0.0750 

23 

21 

0.0 

0.0500 
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APPENDIX B 


In eqn.5.43(a) and 3.43(b) the first term on the R.H.S. 

are equal. The second term on the R.H.S. of 3.43(a) when 

o 

expanded contains terms of the form Sin 6. x.., Sin(6. .-©. .) 

(x^ 1 +Xj^ 2 )j Sin(6 ^^ -f Q ^^ ) x^^-j etc. In order to compare with 

the second term on R.H.S, of 3.43(b) we make use of the following 
inequalities for the defined in eqn. (3.10). 

MJ 



J< 

|Sin 6°^| 


< 

|sui (sy - e^pi (|x.i|+ |Xj,l) 


< 

|Sin ( |xj, 1 


Also we make use of the triangular inequality 

='ir*i2+='i3 i II ^±11 

Finally comparison of like terms in R.H.S. of 3.43 (a) and (b) 


yields 






0 1 o n 

~[ lb. I ICotO^ I |Sin 6 ^ I + I Mj A^lSin( 5 ^r«^^)f]i=j 


k=1 

kji^i 


hi-: 


n^1 


Ai.|Siii(6j^j-e^j) + IIj’ ,1 V'^*3n^l'^5 

Si 


ID 



